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to simulate the rarefied gas flow problems. The DSMC approach is an intuitive and 
probabilistic method which is capable of simulating the low-density and rarefied flows. 
Calculations can be made throughout the rarefied regimes and can overlap with the con- 
tinuum calculations, the degree of overlap depends on the magnitude of the computation 
resources. Because the DSMC computational time is proportional to the number den- 



sity in the computational domain, it is quite expensive to apply the DSMC method for 
the flowfield inside the nozzle. In this research, the DSMC technique is applied to the 
exhausted plume near the nozzle exit and the backflow region. The nozzle exit data, 
required by DSMC, are provided by the nozzle solutions. The nozzle exit flow condi- 
tions are very critical to the accuracy of the plume flowfield prediction since it affects 
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quired about 8 hours of CPU time on the CRAY YMP. It clearly shows that the sonic 
line is shifted to the nozzle exit and substantial portions of the flow field inside the nozzle 
are subsonic. Figure 5 shows the Mach number distribution along the nozzle centerline. 
The comparison of the discharge coefficient between the solution of the Boltzmann-2D 
and the experimental data is shown in Figure 6. The discharge coefficient is defined by 
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Figures 7(a), 7(b) and 7(c) show the Burnett solutions of velocity vectors, local 
Mach number and averaged Knudsen number contours with chamber pressure 0.1 psi, 
under adiabatic wall conditions. In this case, the calculated Reynolds number based on 
the nozzle throat diameter is 40. The averaged Knudsen number with respect to nozzle 



throat diameter, ranges from 0.03 at nozzle inlet to 0.3 at nozzle exit, which indicates 
that the nozzle flow field is in the transitional regime. A decelerating region is found near 
the nozzle centerline from the Mach number contours in Figures 7(b). 

Burnett equations are also solved with the slip-wall boundary conditions. The wall 
temperature is constant. The nozzle chamber pressures are 0.1,0.25,0.5 and l.Opsi. The 
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The calculated pseudo-Reynolds number, 
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10, where the flow inside the nozzle is in the transitional regime. This indicates that the 
Navier-Stokes solutions are not applicable. The merit of adopting the Burnett equations 
as the governing equations is that all the numerical methods for continuum flows are 
applicable. The solution of Burnett equations requires much less computer resources 
compared to the application of DSMC technique. However, the Burnett solutions begin 
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The nozzle chamber pressure is 303 Pa, and the nozzle chamber temperature is 270 K. 
Constant wall temperature was assumed. The inflow boundary is at the nozzle exit. As 
to other boundaries, the vacuum condition is employed, and there is no mass transfer 
across these boundaries. Preliminary results for velocity streamline and density contours 
are plotted in Figures 14. The exhaust plume in the backflow region is clearly visualized 
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Low-Density Nozzle Flow by Solving the Boltzmann Equation”, AIAA 95-0627, 1995. 






Schematic showing of grid setup 
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Figure 3. Mach number contours inside the nozzle for Kn=0.1. 





Mach number contours Inside MSFC A2 nozzle 




Throat Reynolds Number 

Figure 6. Comparison of calculated discharge coefficient with exp 







Velocity Vectors (Burnett) 

Chamber Pressure 0.1 psi, Chamber Temperature 2025 K 



Figure 7(c). Burnctt-2D Knudscn number contours. Adiabatic wall. 



Mach Number at Nozzle Axis 
(Constant Walt Temperature) 



Figure 10. Burnet t-2D AX: Mach number along centerline. 
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(b). Exit Mach Number Distribution 
Figure 12. Burnett-2DAX solution of MSFC-A2 nozzle flowfield. 



DSMC Computational Domain 
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Figure 14. DSMC streamline and Density contour. 


